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ABSTRACT 

Core Accretion (CA), the most widely accepted scenario for planet formation, postu- 
lates existence of '--^ km-sized solid bodies, called planetesimals, arranged in a razor-thin 
disc in the earliest phases of planet formation. These objects coagulate by collisions, 
eventually building planetary cores. In the Tidal Downsizing (TD) hypothesis, an al- 
ternative scenario for formation of planets, grain growth, sedimentation and formation 
of planetary cores occur inside dense and massive gas clumps formed in the outer cold 
disc by gravitational instability. As a clump migrates inward, tidal forces of the star 
remove all or most of the gas from the clump, downsizing it to a planetary mass body. 
Here we consider a rotating and/or strongly convective gas clump. We argue that such 
a clump may form not only the planetary core but also numerous smaller bodies. As 
an example, we consider the simplest case of bodies on circular orbits around the plan- 
etary core in the centre of the gas clump. We find that bodies smaller than ~ 1 km 
suffer a strong enough aerodynamic drag and thus spiral in and accrete onto the solid 
core rapidly. On the contrary, bodies in the planetesimal and larger size range lose 
their centrifugal support very slowly. We consider analytically and numerically the 
fate of these bodies after the host gas clump is disrupted. Planetesimals orbiting the 
protoplanetary core closely remain gravitationally bound to it; these may be relevant 
to formation of satellites of giant planets. Planetesimals on more distant orbits within 
the host clump are unbound from the protoplanet and are set on mildly eccentric he- 
liocentric orbits, generically forming wide rings. These may correspond to debris discs 
around main sequence stars and the Kuiper belt in the Solar System. For the latter in 
particular, TD hypothesis naturally explains the observed sharp outer edge and the 
"mass deficit" of the Kuiper belt. 

Key words: general ~ minor planets, asteroids: general - planets and satellites: 
formation - planets and satellites: methods: numerical - Kuiper belt. 
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1 INTRODUCTION 



1.1 Planet formation: bottom-up or top-down? 



model may be termed a bottom-up scenario for planet for- 
mation. 



An alternative top-down scenario for the origin of plan- 
ets is also physically plausible, although it is currently much 
et al.|[2005l posits that planet formation starts when some less developed (Boley et al. 2010 Nayakshin 2011a|. The 



Core Accretion (CA) theory (Pollack et al. 1996 Alibert 



of the dust grains in the protoplanetary disc grow and sed- 
iment to the disc midplane. This step is then followed by 
a less well understood one in which much larger solid ob- 
jects of a few km or more are made (cf. [Youdin fc Shu|2002[ 



Johansen et al. 20071. These rocky objects, called "plan- 



very first step in this "Tidal Downsizing" (TD) scheme is 
reminiscent of formation of stars, that is the Jeans self- 
gravitation instability of a gas cloud/clump (for example, 
see |Larson||1969| . The cloud is however much less massive, 
e.g., its mass is around the opacity fragmentation limit of 



etesimals", then co-coagulate into terrestrial-like planetary 
cores ( |Safronov||1972[ ). These cores continue to gain mass 
by mergers and accretion of planetesimals. If the cores are 
still embedded in the parent gaseous disc, then a gas at- 
mosphere builds up around the cores. When the solid core 
mass exceeds a few to a few tens Earth masses, the atmo- 
sphere becomes self-gravitating and collapses hydrodynam- 
ically to much higher densities, forming a proto-giant gas 
planet ||Mizuno|[T980} |Stevensonl|1982[ [R"afikov||2006[). This 



lOMj (Rees 1976 Nayakshin 2010b Forgan & Rice 2011 1, 



and is born inside a gas disc orbiting the parent star (e.g.. 
Boss 19971. The second step in TD scenario is similar to 
the first one in CA model, e.g., grains grow and sediment. 
However, this process occurs inside the gas clumps rather 
than the whole of the disc, where the densities are much 
lower (cf . simulations of Cha fc Nayakshin|2011 1 . A massive 
solid core therefore forms inside the gas clump ( |Boss|[T998l 
Helled fc Schubert|[2008l |Nayakshin||2011a[). Rapid inward 
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radial migration of the gas clump (Vorobyov fc Basu]|2005 



20M 



et al. 



Machida et al.||2011| [Michael et al.||2011| |Baruteau| 



2011 1 then expose the clump to the ever increasing 



tidal force of the parent star. Removal of all or a part of the 
original gas cloud by tidal forces of the star may hypotheti- 
cally leave behind both terrestrial-like and gas giant planets 



sequence stars ( Zuckernian|2001 Wyatt|2008 |. Since micro- 
scopic dust in the discs around main sequence dust should 
be blown away rapidly due to radiation pressure from the 
parent stars, a continuous replenishment source for the dust 



is required (Wyatt et al. 20071. The universally accepted 



picture is that the circum-stellar dust results from a frag- 



(Boley et al. 2010 Nayakshin 2010a). Note that while the mentation "cascade" (e.g., Hellyer||1970 Wyatt et aL]|2007 



combination of these physical steps is a very recent devel- 
opment, the TD hypothesis may be viewed as a physics- 
upgraded version of the gravitational disc instability model 
(GI; e.g., |Boss||1997| [Rice et al.||2005| |Rafikov|[200^. The 



complicated fate of gas clumps in a self-gravitating disc was 
also discussed by Mayer et al. (2004!), who also noted that 



these clumps may be tidally disrupted if they migrate in- 
ward. Furthermore, except for the crucial step of the radial 
migrati orP] of the gas clump , the TD scheme is similar to 
ideas oflMcCrea & Williams||l965|. 



The current status of TD planet formation hypothe- 
sis is best described as "work in progress" . Progress is be- 
ing made in terms of building more self-consistent time- 
dependent models of the protoplanetary discs with embed- 
ded massive gas clumps (e.g., [Boley et aL||2011[ [Nayakshin 
fc Lodato|2011 1; addressing the surprising high temperature 
content of Solar System comets ( Vorobyov|2011 Nayakshin 



et al.|p011[ l; chemistry of self-gravitati ng discs (|Ilee et al 
2011[); the "hot" super-Earth planets (|Nayakshin||2011b[); 



and rotation of the Solar System terrestrial planets ( Nayak- 
|shin||2011c[ ). Despite this, no detailed statistical predictions 
in the spirit of population synthesis models of [Ida fc Lin| 
( 20081 |2010l) have yet been made. The TD hypothesis thus 



cannot yet be compared with exoplanetary data in detail. 

One promising way to constrain the TD hypothesis ob- 
servationally is via observations of the early "embedded" 



phase of star (and potentially planet) formation (Dunham 
fc Vorobyov|2012l, especially with the advent of the ALMA 



telescope when resolving individual massive gas clumps in 
large 7? ^ 50 AU discs around young protostars may become 
possible. 



1.2 Can solid debris constrain planet formation 
theories? 

Planets are not the only bodies orbiting their host stars: 
solid bodies from microscopic dust to sub-terrestrial planet 
size objects such as Pluto also form during planet assembly 
both in the Solar System and around other nearby stars. In 
particular, Solar System contains several repositories of ;^ 1 
km sized solid bodies - the asteroid belt, the Kuiper belt, 
the scattered disc and the Oort cloud of long-period comets. 
Starting from the discovery of the infrared excess beyond 
12/im around Vega by IRAS ( Aumann et al.|1984 l, circum- 
stellar dust disc are commonly found around pre- and main 



^ |Donnison &: Williams| (|1975^ showed that the gas clump with 
properties envisaged by |McCrea fc Williania| | |1965[ l , sitting at the 
present location of the Earth, would be tidally disrupted well be- 
fore grains could have grown and sedimented to its centre. Making 
the clump at ~ 100 AU, where the clump may be much cooler, 
and having grains sediment into a massive solid core before bring- 
ing the clump into the inner Solar System by migration is the only 
physically plausible way for the model to work. 



'Wyatt 2008) of larger solid bodies such as comets and aster- 
oids. This view is reinforced by observations of dust resulting 
from such fragmentation cascades in our own Solar System 
( Nesvorny et al.||2003 20101, and also by the signatures of 
coUisional sculpting of the asteroid belt (e.g., Bottke et al 



20051 



It is obvious that the solid debris around stars has a 
rather natural explanation in the context of the CA the- 
ory: these bodies are the planetesimals and "half-grown" 
planetary embryos that did not get consumed or ejected 
completely from the system by the growing planets. The 
question we pose in this paper is this: Can TD hypothesis 
produce solid debris too, and if so, how do the properties of 
that debris differ from the CA model? 

We shall argue below that it is economical and logical 
to have all smaller solid objects to be born inside the same 
"parent" gas clumps where the terrestrial-like solid cores 
grow. We arrive at these conclusions based on the follow- 
ing ideas. In the ID spherically symmetric models of Helled 



and co-authors, and in Nayakshin (|2011a|, only one central 
core is formed. However, ( Nayakshin|201 lb I has shown that 
specific angular momentum of the grains may be too large 
to allow gravitational collapse into just one massive body. 
Gravitational collapse of a rapidly rotating cloud may result 
in formation of not only the central body but also a num- 
ber of smaller objects possibly orbiting the larger one in a 
centrifugally supporting disc. The objects may fragment or 
grow further due to subsequent collisions. In addition, the 
formation of the solid core releases a significant amount of 
energy into the surrounding gas ( |Nayakshin||201l"a[ ), which 
stirs up strong convective motions. This second effect may 
also promote formation of additional solid bodies by gravita- 
tional collapse of smaller grain-dominated regions (cf. §3.6.2 
and 3.6.3 of Nayakshin|2010b I on chaotically oriented orbits. 

The next logical step in our picture is the disruption 
of the parent gas clump, as required by the TD hypothesis. 
We show below that the solid debris population within the 
clump can be divided into two groups: one bound to the 
planetary core (be it terrestrial-like or more massive by that 
point), and the other to the host gas clump. Upon disruption 
of the clump the first group "survives" and becomes planet's 
satellites; the other group gets disrupted and dispersed with 
the gas. Generically, the second population of debris forms 
a "disruption ring" centred on the location of the host gas 
clump disruption. 

Below we study analytically and numerically (for three 
representative cases) the properties of the post-disruption 
orbits of the planetesimal debris. We start in §2 by discussing 
the likely structure of the gas clump before the disruption. 
We show that there are two possible ways of the host clump 
disruption - either tidally or due to an internal energy re- 
lease by the growing protoplanetary core itself. We estimate 
the minimum mass of the core being able to disrupt the host 
clump to be ~ 10 Earth masses. 

In |3.1|we consider the aerodynamic drag acting on the 
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solid bodies within the host clump. We find that bodies 
smaller than ~ 1 cm and larger than ~ 1 km stand a good 
chance of remaining "independent" after the clump disrup- 
tion compared with the bodies between these size limits: 
these suffer strong drag from the gas and must end up join- 
ing the protoplanetary core in the centre of the gas clump. 
In §3.4| we consider analytically the kind of orbits that the 
large solid bodies obtain after the gas clump disruption. 
^ describes the set up of our numerical experiments. Sev- 
eral following sections present the numerical results, with !J5] 
showing the gas flow, §!|6] and IT] focusing on the population 
of unbound and bound solid bodies, respectively, ^contains 
a discussion of the main results of our paper. 



2 HOST CLUMP STRUCTURE AND 
DISRUPTION 

2.1 Setup and terminology 

For definitiveness, we consider a host gas clump of M^c = 
5Mj mass in the numerical part of the paper. We assume 
that the initial density and temperature profiles are that of 
a polytropic sphere with the polytropic index n — 5/2 as 
appropriate for a molecular hydrogen-dominated gas in the 
temperature range from a few hundred K to about 1500 K 
I Boley et al.|2007) ). The clump is initially located at i? = 40 
AU on a circular orbit around the star of mass M« = IMq. 
We shall refer to the host gas clump as such to distin- 
guish it from the solid planetary core, which we frequently 
call "planet". The planet (of mass Mp — 10 M®) is treated 
as a point mass in this paper. We assume that the density 
of the planet is much higher than that of the host clump. 
Therefore, the internal structure of the planet is not affected 
by the host gas clump disruption. We note that this setting 
does not require the planet to consist of high-Z material 
only; volatiles could be present also as long as their densi- 
ties are much higher than the tidal density, M,/2-kR'^. 



2.2 Clump disruption by tidal forces 

We study two limiting cases of clump disruption. In the case 
of disruption due to tidal forces from the star, the host clump 
fills its Roche lobe at the start of the simulation, time t = 0. 
The size of the Roche lobe (Hill's radius) at this location is 



rh 



^r(M^\ 



1/3 



\SMj 



; 4.7 AU 



(1) 



A polytropic cloud with the clump mass Mhc = 5Mj and 
the clump radius, rhc, satisfying rhc = r^ has central tem- 
perature T — 195 K, for reference. 

The mass-radius relation of a polytropic cloud is given 
by 



The oc M^-" oc i\4/ for n = 5/2 



(2) 



This implies that the host clump expands rapidly as mass is 
lost. We therefore expect a prompt tidal disruption of such 
a host clump by Roche lobe overfiow: The increases while rh 
decreases as Mhc drops. 

We shall add to this that the stabilising effect of the 
host clump outward migration due to the mass exchange 



between the gas clump and the star, discovered by |Nayak^ 
shin & Lodato (20111, does not occur in our simulations 



because the clump is destroyed rather rapidly in all of the 
simulations below. In addition, we find that disruption pro- 
ceeds via both LI and L2 points, in contrast to what is found 
for much more gentle "hot disruptions" at 7? ~ 0.1 AU sep- 



arations by Nayakshin & Lodato (20111. For both of these 



reasons there is no outward torque on the host clump; the 
tidal disruption has a runaway character. 

2.3 Disruption by an internal energy release 

In the second class of models we consider here, the clump's 
initial radius is smaller than the Hill's radius Vh- In par- 
ticular, for the two simulations presented below, we set the 
initial central temperature to T = 500 K, which corresponds 
to the clump's initial radius of rhc = 1-84 AU. We then as- 
sume that a sudden burst of energy release occurs in the 
centre of the clump. Physically, we relate the burst to the 
assembly of the planetary core, as is found in ID simula- 
tions of |Nayakshin| ( [2011a[ ). In detail, the energy released by 
the solid core is passed to the surrounding gas by radiative 
diffusion. If the radiative diffusion time scale is longer than 
the planetary core's assembly time and the injection energy 
is large enough then even an isolated gas clump may be dis- 
rupted. An example of this is simulation M0a3, Figures 5 to 
7 in [Nayakshin] ( |2011a| ). 

We estimate the energy released by the planetary core 
to be about its binding energy. 



2r„ 



= 5 X 10*" erg (' ^'^^ 



V lOA/ff 



5/3 



P^ 



(3) 



where G is the gravitational constant, and pp is the solid 
core density in g cm~^. 

Let us now compare the core's binding energy with 
the minimum amount of energy needed to disrupt the host 
clump. The total energy of a polytropic gas clump, She is 
given by 



-Ehr, ~ — 



3 - n GMl^ 
5 — n The 



(4) 



For simplicity we assume that the clump remains polytropic 
after the energy injection (which is not obvious at all; we 
will come back to this later) with same n but a different 
adiabatic constant. The energy input sufficient to disrupt 
the clump, Af/o is the energy needed to inflate the clump to 
the point of its tidal disruption, i.e., when rhc increases and 
becomes equal to the Hill's radius, r^. Thus, 



5 — n V Th 



GM, 



^hc 



(5) 



For reference, GM^^/rhc ~ 2.4 x 10'*^erg for clump mass 
Mhc = 5Mj and clump radius rhc ~ 1.84 AU. 

As energy released by the core heats the surrounding 
gas, it is instructive to compare the required energy injection 
AUo with the initial total internal energy of the host clump, 
Uo: 



AUo ^ _ /, 3 

Uo 



1- 



nic 

Th 



(6) 



In the limit of the initially small host clump, rhc ^ '"h, the 
disruption energy is AUo ~ 0-2Uq for n = 5/2. In the case 
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under consideration, Hic = 1-84 AU, and r^ = 4.7 AU, and 
thus At/o = 0.123C/0. 

From the foUowing we conclude that the energy required 
to unbind the host gas clump is of the order of 3 x lO'*" erg. 
This is of the same order as the binding energy of the plane- 
tary core of mass Mp = lOM^. Therefore, we see that, to dis- 
rupt the host clump with the energy released by the plane- 
tary core, the mass of the core must be at least Mp ~ 10 Af^. 
We also note in passing that the similarity of this mass to 
the critical solid core mass in the CA theory (e.g., [M"izuno| 
1980 Stevenson|1982 \ seems to be a pure coincidence. While 
the minimum mass capable of disrupting the gas host clump 
is the function of the conditions in that clump (mass, age, 
opacity, etc.), the CA critical mass is a function of the core's 
surroundings - the location from the star, opacity and the 
planetesimal accretion rate from the surrounding protoplan- 
etary disc (e.g., Rafikov|[2011 1. 



3 SOLIDS WITHIN THE CLUMP 



Nayakshin (2010b 201 la I finds that when grains contained 



in the host clump grow by the hit-and-stick mechanism to 
sizes of the order of s ~ 10 cm, their sedimentation within 



~ 1000 years becomes possible (see also McCrea & Williams 
|1965[ [Boss|1998| . Accumulation of these to the centre of the 
clump creates a region dominated by grains rather than by 



gas. This region, called "grain cluster" by Nayakshin ( 2010b I 



then becomes gravitationally unstable and collapses to form 
a massive core composed of high-Z materials. 

As argued in §1.2[ due to rotation and chaotic convec- 
tive gas motions in the centre of the host clump, one may 
expect numerous smaller bodies to form in the centre of the 
clump as well. It is the fate of these bodies that interests us 
here. 



3.1 Aerodynamic drag in the host clump 

We shall now consider what happens to solids of different 
sizes, s, from microscopic grains to asteroid-sized bodies, 
if they are placed within the host clump. As we shall see, 
due to the aerodynamic drag that these solids suffer, small 
grains (s ^ 0.1 — 1 cm) are nearly frozen in with the gas, so 
have to follow its motion, whereas bodies larger than s ~ 1 
km experience negligible friction with the gas. Bodies in the 
intermediate size range are most likely to end up joining the 
solid core. The smaller and the larger objects may survive 
the gas clump disruption and become either planet satellites 
or independent bodies in heliocentric orbits. 

To quantify this discussion, we model the host gas 
clump in this section only as in the analytical model of 
|Nayakshin| ( |2010a[ ), which shows that the outer radius of 
the host clump is approximately independent of the clump's 
mass, The = 0.8 AU kj (10*yr/thc)^ , where ihc is the 
host clump's age, and fc* is dimensionless opacity (cf. §2.1 
of Nayakshin||2010a I . The age of the host clump chosen for 
the representative calculation below is 10* years, the mass is 
lOMj, and opacity A;» = 1. Note that in detail the structure 
of the host clump is different from that of the gas clumps 
we consider in the numerical of the paper, but the main 
conclusions that we reach in this section are qualitatively 
unaffected by this. 



For simplicity, clump rotation is neglected in this sec- 
tion, but we note that results for rotating gas clumps are 
actually quite similar as long as gas pressure is the main 
means of support against gravity for the host clump (in the 
opposite case the gas clump would be unstable to various 
fluid instabilities; see, e.g., chapter 7 of [Shapiro fc Teukol-| 
[sk^p8 3). 

Our goal here is to calculate how long it takes for a 
grain of size s and a given initial condition for the grain 
position and velocity to settle to the centre of the clump due 
to aerodynamic friction with the gas. To achieve this goal, 
we solve for the grain's motion within the host clump using 
the formulae of |Weidenschilling| ( |1977[ ) for the aerodynamic 
forces acting on the grains. This allows us to calculate the 
radial sedimentation velocity of the grain, Vr, and deflne the 
sedimentation time as 



r 

\Vr\ 



ij) 



where r is distance to the clump's centre. For solids released 
from rest results of such calculations are presented in fig. 2 



of Nayakshin et al. (20111. Here we present a very similar 



calculation for dust grains that are on initially circular orbits 
around the host gas clump's centre. 

Figure IT] shows the sedimentation time scale for grains 
released on circular orbits with initial radii r of 0.02, 0.1 
and 0.5 times rhc for the blue dot-dashed, brown dashed 
and black solid curves, respectively. For comparison, the red 
dash-triple-dot curve shows the grain sedimentation time for 
a grain released from rest, as in Nayakshin et al. (2011 1. 



Without rotation, grains larger than a few cm manage 
to sediment to the centre of the gas clump within its age. 
Smaller grains are tightly bound to the gas, and if the host is 
disrupted, these grains are released into the disc around the 
protoplanet. The grains closest to the solid core could have 
been thermally reprocessed into crystalline materials; mix- 
ing these with surrounding disc's ices and incorporating into 



comets may explain their puzzling compositions (Vorobyov 
|2011[ [Nayakshin et al.|2011| ). 

In a purely spherical geometry with no rotational sup- 
port for grains, large grains (s ^ 1 m) fall into the centre 
and join the solid core on the free-fall time, which for a con- 
stant density model clump is constant with radius (cf. the 
horizontal part of the red dash-triple-dot curve) . 

For grains on circular orbits, however, centrifugal force 
balances gravity, and the grains sediment only because of a 
gradual angular momentum loss due to aerodynamic forces. 
These forces become progressively less important for larger 
grains, and therefore the grain sedimentation time increases 
with s again (cf. the black, the brown and the blue curves in 
Fig.[I|. As a result, grains larger than 1-10 km may orbit the 
centre of the clump for times comparable with the clump's 
lifetime (which we assume of the same order as the clump's 
age here, as presumably the clump continues to migrate in- 
ward at roughly same migration speed). 

Concluding, we see that solid objects larger than a few 
km in size do not necessarily contribute to growth of the 
solid planetary core as it may take too long for these bod- 
ies to spiral into the core: the host clump is likely to be 
disrupted before such an inspiral occurs. 
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Figure 1. Sedimentation times versus grain size for grains set 
initially on circular orbits, except for the red dash-triple-dot curve 
which is for a grain released from rest. 



3.2 On planetesimal birth and sizes 

In this paper we do not simulate the phase of planetesi- 
mal formation due to numerical challenges. Previous simu- 
lations of the whole proto-planetary disc by |Cha fc Nayak-| 
shin (20111 lacked numerical resolution (to resolve smaller 
objects) and resulted in formation of single self-bound mas- 
sive (M ^ a few M^) clusters of dust particles inside the 
gas clumps. These were held from a further self-gravitational 
collapse artificially by employing a finite gravitational soft- 
ening length in the simulations. Increasing numerical resolu- 
tion further (that is, the number of SPH and dust particles) 
by up to several orders of magnitude is needed to resolve 



small, e.g., ^ 0.01 AU scales, on which planetesimals we 
are interested here could form. This is beyond our present 
numerical capabilities. 

Figure [T] shows a difficulty for formation of planetes- 
imals in the TD scheme, very much similar in nature to 
the well known "metre-size barrier" for the planetesimal 



growth in the protoplanetary discs ( Weidenschilling,, 1977) . 
In the CA theory, aerodynamic coupling with the gas makes 
~ metre-sized boulders migrate radially inward in the disc, 
so that they are probably lost to the protostar before they 
become planetesimals. 

Similarly, Figure [l] shows that solids of intermediate 
sizes migrate radially inward very rapidly even if they are ini- 
tially set on circular orbits around the centre of the host gas 
clump. Such objects thus must join the massive solid core 
there. To arrive at km-sized and larger bodies that suffer 
much weaker aerodynamic drag, smaller grains must thus 
somehow "jump" from being small to being large without 
suffering the aerodynamic drag. 

We leave a detailed study of this issue to a future paper, 
but make some suggestions on how this difficulty may be 
avoided. First of all, again reminiscent of the well known 
ideas in the CA theory context, channels for a rapid growth 
of solids may be available: via self-gravitational instability 



( Safronov|19"72l|Goldreich fc Ward|1973l l of dust-dominated 
regions (although not in necessarily in the shape of a disc) , 
or via turbulence (e.g., Youdin fc Goodman|2005 Johansen 



et al.|2007||Cuzzi et al.[2008 i' 



Secondly, and unlike the CA "metre-size barrier" prob- 
lem, here it is much less clear that aerodynamic drag does 
make all the medium-sized solids to accrete onto the central 
object. In the case of a proto-planetary disc (the CA model), 
the protostar contains almost all the mass and is thus the 
natural central (nearly) point-mass around which the disc 
rotates. Grains migrating inward must therefore end up in 
the star. In the TD case, however, the mass of the "mas- 
sive" solid core we are thinking about is measured in Earth 
masses, which is ~ 10"'' of the total host gas clump mass. 
If there are turbulent or convective motions in the inner re- 
gion of the clump (which is almost a certainty), then the 
solid core itself will be affected by the gas motions through 
its coupling to the gas via gravity of the latter. Thus it may 
not "sit" in the centre of the clump accreting all the medium- 
sized solids sedimenting there. Furthermore, motion of the 
grains would also be affected, with turbulent and convective 
motions of gas dragging the grains around. There is clearly a 
limit here; too much of turbulence would prevent grains from 
sedimenting into the centre in the first place, but moderate 
turbulent motions may be expected to delay grain accretion 
onto the solid core and lead to formation of new condensa- 
tion centres which perhaps would lead to formation of the 
larger planetesimals that we study in the rest of the paper. 



It would also be desirable to understand just how large 
the planetesimals formed inside the host gas clump are likely 
to be. Following §3.6.2 and 3.6.3 of [Nayakshiiil ( |2010b[ ), 
we find that the linear extent of the region of the gas 



cloud within which the first gravitational collapse of grain- 
dominated regions could occur is about ~ 0.1 times the size 
of the gas clump. For the total grain mass in the collapsing 
region of 10 M^ , the expected velocity dispersion of frag- 
ments formed by collapse is then ~ 0.5 km s^. 



Now, collisions of large planetesimals of equal size, a, 
would lead to their fragmentation if the collision velocity (~ 
the velocity dispersion calculated above) is larger than the 
escape velocity from the surface of the planetesimal (e.g., 
Stewart & Leinhardt 20091. This would then suggest that 



collisions of equal size bodies inside the "grain sphere" would 
split planetesimals smaller than about 1000 km, while larger 
bodies would "stick" . One however need to estimate the fre- 
quency of such collisions. Assuming that the grain sphere 
fragments into bodies of approximately equal size, we find 
that all the planetesimals would have had at least one catas- 
trophic (shattering) collision within a thousand years inside 
the grain sphere for a ;$, 100 km. Collisions with smaller bod- 
ies - fragments of the "original" planetesimals - increase 
this collision rate estimate further. We thus conclude that 
expected sizes of planetesimals surviving the dense environ- 
ment of the inner region of the clump are at least a few 
hundred km. Presumably both smaller and larger solid ob- 
jects could then be obtained over much longer time scales 
by coUisional evolution of the disrupted population. 
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3.3 Planet's influence radius: bound and unbound 
debris 

Having understood the range of body sizes (s ^ Ikm) that 
may be considered as independent for the duration of the 
host clump, we now shift out attention to what may happen 
when the host clump is disrupted. In detail the answer on 
this question is quite complicated as it depends on the na- 
ture of planetesimal orbits within the clump and the exact 
way the clump is disrupted. However, for not too eccentric 
orbits, we can divide the planetesimal population on those 
that are "close" and those that are "far" from the solid plan- 
etary core. The former have a fair chance of remaining bound 
to the solid core when the host clump is disrupted, and the 
latter become unbound independent objects. 

To quantify this, consider the simplest setting, placing 
a solid core of mass Mp in the centre of the gas clump and 
setting planetesimals on circular orbits around the planet. 
The rotation of the planetesimal's disc is assumed prograde 
with respect to the orbit of the host clump around the star. 
The circular speed of the planetesimals around the planet is 
given by 



G 



Mp + M{r) 



(8) 



where M{r) is the gas mass enclosed inside radius r (distance 
from the centre of the host clump) . 

If planet's mass is low enough to build up a massive self- 
gravitating atmosphere around it (which is expected to be 
the case for an adiabatic young gas clump where the critical 
core mass is above 100 Mq ; see Perri & Cameron 1974 1, then 



the gas density is nearly constant in the clump's centre. Let 
us call that central density po; thus M{r) « (47r/3)/9or^. 
We can now define the planet influence radius, Vi, such that 
M{ri) = Mp. Apparently, 

1/3 

(9) 



3Mp 



47rpo_ 

Note that since r-hc = (3Mhc/47rpmcan)^^^, where pmoan is 
the mean density of the host clump, then we have 



' r-h 



( 



M„\i/3 



\Mhc 



po 



1/3 



0.18 riK 



po 



1/3 



(10) 



Solids inside r^ are gravitationally bound mainly to the 
planet, whereas those outside are bound by the enclosed gas. 
We expect that when the host clump is disrupted, objects 
within ri remain bound to the solid planet whereas objects 
outside this radius become unbound from the planet. 

3.4 Post-disruption orbits for unbound debris 

We now build an approximate analytical theory for the or- 
bits of the planetesimals after the host clump tidal disrup- 
tion. To this end we assume that gas clump disruption is 
instantaneous. We also assume that solids within r^ remain 
bound to the planet and continue to follow its motion around 
the central star, and we concentrate here only on solids out- 
side the influence radius. 

Before the disruption, the position of a planetesimal 
with respect to the star is given by R -I- r, where R is the 
vector connecting the star with the centre of the gas clump 
(solid planet's location), and r is the vector connecting the 



centre of the host clump with the given planetesimal. The 
instantaneous velocity of the planetesimal is also comprised 
of two parts, one due to the heliocentric Keplerian motion 
of the gas clump, V, and the other given by the circular 
prograde rotation velocity, iidrc, around the centre of the 
host clump. 

We now call V + Av the planetesimal's velocity right 
after the host clump disruption. Physically, Av is the ve- 
locity with which planetesimal becomes unbound, and may 
be called a velocity "kick". We expect that < Av < Ucirc 
for the energy conservation reasons. 

We shall now calculate the new heliocentric orbit of 
the planetesimal by considering its specific energy and an- 
gular momentum. The specific energy and specific angular 
momentum of the host gas clump before the disruption are 
Eo = — GM«/27? and Lo = R x V, respectively. The post- 
disruption specific energy, E, and specific angular momen- 
tum, L, of the planetesimal are 

E = -^M^ + 1|V + Av|^ (11) 



|R + r| 2 
(R + r) X (V + Av) 



(12) 



where, M« is the mass of the central protostar. Decomposing 
these expressions into Taylor series with respect to small 
parameters r/R and Av/V , up to the linear terms only, and 
writing E « Eq + SE, L^ « {Lo + 5LY (the latter is possible 
since Av is in the clump's orbital plane by setup), we have 



5E = V 



i — 



r V- Av 



(13) 



From this we see that the semi-major axis of the planetesi- 
mal's orbit after disruption is 



GM, 



R\ 1 + 



5E 

~E^\ 



R 



l-h2 



Similarly, we can write 

SL ^ rxV + Rx Av . 

Now, using vector identity 

[A X B] • [C X D] = (A • C)(B ■ D) 

we find that 

LoSL =1/^R • r + R^V ■ Av 



K 


R r 


V- AvXl 
V2 )\ 


R^ ' 






(15) 


(A 


D)(B 


C) , (16) 



(17) 



Comparing this expression with equation|13| we observe that 



5E 

Eo 



,5L-Lo 



^0 






(18) 



where we also recalled that the orbital plane of the planetes- 
imals does not change since the proto-disc of planetesimals 
rotates in the prograde direction in the orbital plane of the 
host clump by assumption. 

The eccentricity of the planetesimal's orbit is 

2EL^ 



"" ^ ^ G^Mi ' 

which can be decomposed to become 



SE 5L SE5L 
Eo Lo Eo Lo 



(£)' 



(19) 



(20) 



Using equation [18] we find that the linear terms cancel, and 
the always positive result is 



.(14) 
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(21) 



also vary, likely increase, if planetesimal orbits in the pre- 
disruption clump are eccentric. 



3.4- i Eccentricity - semi-major axis correlation 

We can eliminate 5E and Eq in favour of orbital elements. 
Since Eo = -GMJ2R, and a = Glvh/2\E\, equation [2l] 
actually shows that 



sv^'' M 



(i) 



ao 



(22) 



where ao = R, and 5a ^ a — R is the difference in the semi- 
major axis of the planetesimal's post-disruption orbit and 
that of the host gas clump before disruption. 

Equation [22] predicts a correlation in the orbital ele- 
ments of planetesimals after disruption: the belt of planetes- 
imal debris left after the host gas clump disruption should 
have circular orbits at the centre, a = ao, and increasingly 
eccentric orbits towards the belt's edges. 

We note one rather interesting feature of equation |22| it 
does not depend on the properties of the host gas clump or 
the solid protoplanetary core inside. The eccentricity - semi- 
major axis correlation should thus be a general property of 
the debris rings left after disruption of the host gas clumps 
(as long as the pre-disruption orbit of the clump is nearly 
circular). 

We also note that for the case of the proto-disc of plan- 
etesimals not coinciding with the orbital plane, the inclina- 
tion of the planetesimal orbits, i, after the disruption is non 
zero but can be shown to be small. 



3.4-2 Rings versus discs 

Despite the universal correlation we just discussed, the prop- 
erties of the host gas clump are still imprinted on the range of 
possible planetesimal orbits through the fact that the pre- 
disruption protoplanet is finite in spatial extent, and that 
the maximum velocity kick to a planetesimal released by 
the clump disruption, A«, is also finite. Referring to equa- 
tion [14] we see that the width of the "disruption ring" - the 
ring occupied by the planetesimals left over after the clump 
disruption - is 



4 max 



R r V Av 

+ 



(23) 



ao R^ V^ 

Noting that Au ~ Hdrc oc r within a constant density host 
clump, and that the maximum possible value of r is the 
Hill's radius of the host clump, rn, we can summarise this 
prediction by writing 



4C- 



ao 



(24) 



where (^ is a dimensionless number probably smaller than 
but comparable to unity. For the setup of this paper, in 
particular, where ru ~ 4.7 AU and ao = 40 AU, we have 



0.5C 



(25) 



We shall see below that ( ~ 0.7— 1 for the three simulations 
performed in the paper. In general this parameter will de- 
pend on how concentrated the distribution of planetesimals 
is within the host protoplanet before the disruption, with 
more compact distributions leading to smaller ^. It would 



4 NUMERICS 



4.1 Method 



We now turn to hydro/N-body simulations for a more de- 



tailed investigation of the problem. SPH ( Gingold & Mon- 
laghan 1977 Lucy 1977[ ) is a Lagrangian simulation algo- 
rithm well suited for irregular and self-gravitating systems, 
SPH has been applied to a variety of astrophysical con- 



texts (Monaghan 1992 Springel 20101. In this paper, we 



use GADGET-3, an updated version of the SPH/N-body code 



presented by Springel ( 2005 1 



In both of the disruption scenarios that we study in 
this paper, host clump disruption occurs on a time scale 
shorter than the Kelvin-Helmholtz time of the clump. An 
adiabatic equation of state for the gas is thus sufficient for 
our purposes. As hydrogen is molecular for temperatures 
smaller than ~ 2000 K, we chose the polytropic index, n, 
to he n — 5/2, which corresponds to the ratio of specific 
heats of 7 = 1.4. The number of SPH particles in each of 
the simulations presented is Ngph = 10®, so that the mass 
of each particle is 5 x W^'Mj. 

The protostar, the planet and the planetesimals are all 
modelled as N-body particles of appropriate masses. Accre- 
tion onto the planet is not allowed, but accretion onto the 
protostar is simulated with the sink particle approach (as 
in |Cha fc Nayakshin|20lT] ). This is necessary to prevent the 
build up of very short time step SPH particles near the star 
due to a non-negligible artificial viscosity of the code in this 
low density region, which is additionally of little interest to 
our study. 

The N-body particles interact with gas only through 
gravity (except for gas accretion as described above). We 
set the total mass of the planetesimal disc to Md = O.IM^. 
This mass is shared equally between JVpi = 2 x lO** particles. 
The mass of each is thus ~ 3 x lO'^^ g, and the corresponding 
radius is about 200 km. For such low masses, gravitational 
interactions between planetesimal particles and their effect 
onto the planet and the gas are negligible, even though these 
interactions are included in the calculation. The large linear 
size of the planetesimals allows us to neglect the aerody- 
namic coupling between them and the gas completely. 



4.2 Initial conditions and disruption method 

Below we present three numerical simulations. The first of 
these, labelled UO, uses the initial condition for the poly- 
tropic gas cloud described in §2.2| In this case the host gas 
clump exactly fills its Roche lobe at the beginning of the 
simulation. The clump is thus disrupted due to tidal forces 
of the star only. 

In the other two cases, the host clump radius, rhc ~ 1.84 
AU, is initially smaller than the Hills radius, r^ ~ 4.7 AU. 
The host clump central temperature is 500 K for this initial 
condition. As explained in §2.31 we assume that at i = the 
central solid core releases a given amount of energy which is 
then instantaneously added to the internal energy of the gas 
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particles within radius rcj = 0.2 AU (chosen somewhat ar- 
bitrarily). The inner region is then over-pressured compared 
to its surroundings, which drives an expansion that puffs up 
the whole gas clump. 

In §2.3| we found that injecting energy A[/o = 0.123[/o) 
where Uo is the total internal energy of the gas protoplanet, 
should increase rhc to rh- This derivation assumed that the 
clump keeps its polytropic structure even after the energy 
injection, which is clearly not actually correct for the equa- 
tion of state we use. We therefore performed two simulations 
with energy injection bracketing the critical injection energy 
AC/o. These are labelled UIO and U15, so that AU = 10% 
and %15 of Uo for the two simulations, respectively. We find 
that both of these simulations resulted in the disruption of 
the host gas clump, although the process was much faster in 
U15 than in UIO. There are interesting differences between 
the orbital parameters of the planetesimals in UIO and U15, 
making both simulations worth presenting here. 

The exact outcome of the gas protoplanet disruption de- 
pends on a number of assumptions about the pre-disruption 
stage, such as (i) the exact location of the planet at the 
moment of the gas protoplanet disruption (which needs not 
be the centre of the gas clump in general); (ii) the duration 
of the disruption process (rapid or slow); (iii) the distribu- 
tions and orbits of solids within the clumps, and their total 
mass, which may start influence the dynamics of planetesi- 
mals if their total mass is comparable or larger than that of 
the planet, and (iv) the orbit of the gas clump around the 
protostar. 

The parameter space of the problem is too large to cover 
in this first study. Therefore, we only aim here at learning 
the most roburst features of the results, which we hope will 
be generally correct within a factor of a few despite all the 
uncertainties pointed out above. 

The planetesimal particles are placed in a geometri- 
cally thin disc on circular orbits around the centre of the 
gas clump with prograde velocity given by equation [8] The 
outer radius of the planetesimal disc, r^, is 0.4 AU for the 
simulations U10,U15 and 1 AU for simulation UO, respec- 
tively (to scale approximately as 0.2rhc for all the cases). 
The inner radius of the planetesimal disc, r^, is set to 10% 
of rd- Introduction of the inner disc cutoff at rd does not 
change our results at all but speeds up simulations consid- 
erably. Particles at r ^ Tin are bound very strongly to the 
solid core (planet), so that disruption of the host gas clump 
hardly changes their orbits. These "uninteresting" planetes- 
imals have short time steps and are expensive to simulate. 

The planetesimal disc has the intitial disc surface den- 
sity profile of Spi(»') oc r~^ . As partciles in the disc do not 
interact significantly due to their low masses, behaving as 
test particles, such a choice for Spi(r) has no consequences 
for the results but allows an approximately uniform param- 
eter sampling in terms of the initial distance r to the solid 
planet. 

Initially, the gas host clump rot ates arou nd the pro- 
tostar with the Keplerian velocity, yjGMt,/R. The orbital 
period of the initial orbit, Po, is ~ 253 yrs. We use Po as a 
unit of time in presenting the results below. All of the sim- 
ulations were run until time t — 10, i.e., 10 initial orbits of 
the clump. 



5 DYNAMICS DURING HOST 
PROTOPLANET DISRUPTION 

5.1 Simulation U15 

In this simulation, the energy injected into the gas exceeds 
the critical energy needed to disrupt the host clump (cf. 
eq. |6|. Therefore, a rapid disintegration of the gas clump 
is expected. Fig. [2] shows the projected gas column density 
at two early stages of the simulation, t = 0.06 (left panel) 
and t — 0.25 (right panel). The projection is done along the 
z-axis, e.g., along the direction normal to the host clump's 
orbital plane. The figures are centered on the position of the 
solid core (planet), shown with the thick green symbol at 
the centre of each panel. Vectors show the gas velocity field, 
likelywise centered on the solid planet's velocity. The black 
dots are individual planetesimals. 

Initially (left panel), clump expansion appears nearly 
spherically symmetric even though the outer gas shells did 
overfiow the Roche lobe by that time. However, the flow be- 
yond rh is controlled more and more by the star's tidal fleld. 
This is evident in the right panel, where the flow becomes 
highly assymetric at large r. The star is positioned exactly 
North of the clump in the right panel. One observes a flow 
of gas towards the star at the upper left corner of the panel, 
and away from the star in the lower part of the panel. The 
planetesimal disc is also becoming affected, visually follow- 
ing the motion pattern of the gas. 

The host clump disruption is indeed dynamic as the 
right panel corresponds to time just slightly later than one 
dynamical time, defined as Po/2n. 

Fig. [3] shows the same as Fig. [2] but now at times t = 1 
and 10 in the left and the right panels, respectively. The fig- 
ure presents a larger field of view centered onto the star now. 
The velocity field is also centered onto the star rather than 
the planet. Note in the left panel that the planetesimals dis- 
rupted off the initial gas clump continue to follow orbits sim- 
ilar to that of the densest regions of gas. This is hardly sur- 
prising given that potential due to the host clump initially 
dominates the planetesimals' orbits (until the host clump is 
completely disrupted). At late times (see the right panel), 
the host clump is disrupted into a ring centered about the 
initial clump's separation from the star. The planetesimals 
are also spread into a ring-like feature which shows several 
streams. The streams appear to have some eccentricity to 
them. We also note that the fine structure of the streams is 
now not correlated to the gas component. This is probably 
caused by the gas gravitational potential being smoothed 
out due to circularisation of the gas flow, leading to a re- 
duced gas gravitational force on the planetesimals. 



5.2 Simulation UIO 



As derived in |2.3[ the critical injection energy sufficient 
to inflate the clump to the point of its tidal disruption is 
AUo = 0.123f/o for the parameters of the clumps we con- 
sider. Therefore, if the approximate theory of §2.3| were cor- 
rect, the clump would not be disrupted at all in this par- 
ticular simulation as AU = O.lUo < AUo- However, the 
argument given by equation [6] assumes that the energy in- 
jected into the central regions of the gas is shared by the 
whole clump, so that gas in the clump effectively finds it- 
self on a new (single) polytropic relation. In the simulations. 
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Figure 2. The top view of the gas surface density map for simulation U15 at time t=0.06 (left panel) and t = 0.25 (right panel). The 
map is centered on the position of the solid planet (core), which is shown with the thick green symbol. Vectors show the velocity field, 
also centered on the velocity of the planet. The black dots show locations of individual planetesimals. 
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Figure 3. Same as Fig. [2]but at later time; i = 1 (left panel), and i = 10 (right panel), except that the center of the coordinate system 
is now the protostar (the red symbol at the center of the panels). Note that planetesimals initially continue to follow the orbital motion 
of the densest parts of the gas streams at t = 1 but separate out at t = 10 when the gas streams are dispersed and circularised into a 
ring. 
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R [AU] 

Figure 4. The internal density structure of tlie expanding pro- 
toplanetary clump for simulation UIO {AU/Uo = 0.1) at time 
t = 0.05. The red solid line is the actual gas density from the 
simulation, while the black solid line is a polytrope of 5Mj mass 
and 4.7AU radius. The latter curve is expected based on the ap- 
proximate theory shown in j ]2.3| The figure shows that a single 
polytropic relation approximation breaks down after the energy 
injection, although the approximation is still useful in roughly 
determining the critical injection energy AUo- 



however, the central regions are shifted to a different poly- 
tropic relation by the energy injection, whereas the rest of 
the clump remains at the initial one. One may thus expect 
the simulation results to differ from the simple analytical 
prediction. 

Figure HI shows the radial density profile of the host 
clump (red curve) in simulation UIO at time i = 0.05 and 
the theoretical polytropic density profile (black curve) with 
n = 5/2 inflated by the energy input from the solid core as 
assumed in §2.3[ One notices significant differences between 
the two curves, with a tail of gas material extending further 
out in the red curve. It is this tail that allows some material 
to siphon out of the Roche lobe and for the gas clump's 
radius to continue swelling (cf. equation [2] on this point) as 
mass is lost. 

The left panel of Figure [5] illustrates this initially slow 
expansion. The host clump is largely intact at i = 3 in UIO. 
This is in stark contrast to simulation U15 in which the host 
gas clump was completely obliterated by t = 1 already (left 
panel of Fig. [3|. However, the host gas clump is eventually 
destroyed in simulation UIO also. The gas column density 
and the distribution of planetesimals appear to be quite sim- 
ilar at late times in UIO and U15 (t = 10, right panels of 
Figs. [5] and [3] respectively). We shall make a more detailed 
analysis of planetesimals' orbits in !|6] below. 



5.3 Simualtion UO 

In this simulation, the gas protoplanet is more extended 
initially, so that rhc = r^ at i = 0. No energy input from 
the solid core is assumed. Figure [6] shows two snapshots for 
simulation UO at times i = 4 (left panel) and i = 10 (right 
panel) in the same format as Fig. 15] The disruption process 
is not as rapid as in simulation U15 but is a little faster 
than in UIO. The morphology of the gas flow is different: 
the disrupted gas spiral at t = 4 is wider in UO than it is in 
Fig-i 

Since the initial disc or planetesimals inside the gas pro- 
toplanet is large to begin with (vd = 1 AU for UO) than for 
simulations UIO and U15 (ra ~ 0.4 AU in both), the flow 
of disrupted planetesimals is initially wider in UO (see the 
banana-shaped feature in the left panel of Fig. [6]) The end 
result is at least visually not too dissimilar from runs UIO 
and U15, however. 



6 ORBITS OF THE UNBOUND 
PLANETESIMALS 

We now turn to analysis of the unbound population of plan- 
etesimals in terms of their orbital elements at the end of the 
simulations. To avoid confusion, tt must be stressed that, 
all of our planetesimals remain bound to the central star 
at the end of the simulations, but some are also bound to 
the planet, as planetary satellites. The "unbound" popula- 
tion of planetesimals are those particles moving on their own 
independent heliocentric orbits. 

In the simulations, the bound and unbound populations 
are differentiated according to two conditions. We deflne the 
specific energy of a planetesimal with respect to the planet. 



Er, 



GM„ 



+ 1 (v-Vp 



(26) 



where r is the distance between the planet and the planetesi- 
mal, V and Vp are the planetesimal and the planet velocities, 
respectively (note that these quantities are calculated here 
at i = 10 whereas in §2.3| we measured relative positions 
and velocities of planetesimals and the planet at i = 0). The 
planetesimal is considered unbound if 



Ercl > , 

and bound if 

Erci < and r ^ rf^ 



(27) 



(28) 

where r'^ = i?o(Mp/3M*)^/^ is the HiU's radius of the solid 
planet. 

The two populations should be analysed differently, of 
course. The orbits of the unbound population should be de- 
fined with respect to the central star, whereas the orbits of 
the bound planetesimals are best defined with respect to the 
planet. In the rest of this section we consider the unbound 
part of the planetesimals only. 



6.1 Simulation U15 

Fig. JTI shows the orbital eccentricity (bottom panel) and the 
inclination (top panel) of the unbound planetesimals versus 
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Figure 5. Gas column density projections similar to Fig. H] but for run UIO at time t = 3 (left panel) and t = 10 (right panel). As the 
planet is less extended than in test U15, the disruption process is much more gradual. At t = 3 (3 whole orbits of the host clump around 
the star), only a small fraction of gas mass is lost through the Roche lobe overflow. No planetesimals have yet been unbound from the 
clump. However, at late time (the right panel), the structure of the disrupted gas and planetesimal population is quite similar to that in 
U15. 
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Figure 6. Same as Fig. [s] but for simulation UO, the case where the initial cloud is cooler and already fills its Roche lobe. The 
corresponding times are t = A (left panel) and t = 10 (right panel). 
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the semi-major axis of their orbits (a) at t = 10. As ex- 
plained above, their orbital parameters are calculated with 
respect to the star as their orbits are hehocentric. 

The top panel of figure [T] shows that planetesimals con- 
tinue to have very small orbital inclinations after the disrup- 
tion, with the mean inclination angle of only i « 1°. This 
is natural since in our setup tidal disruption of the host gas 
protoplanet is symmetric around its orbital plane. In fact, 
the approximate analytical theory of §2. 3| predicts that in- 
clination of planetesimal orbits should remain exactly zero. 
Therefore, the inclination found in the simulation must re- 
sult either from a small but finite asymmetry with respect to 
z -^ —z inversion developing during the disruption, or due 
to effects not taken into account in the analytical theory. 

The central a ~ 32 to a ~ 40 AU region of this plot 
shows a higher dispersion in the values of i than do the 
more distant regions. To understand the origin of this result, 
we note that the central region is the one where most of the 
gas ends up after being disrupted, e.g., see the right panel of 
figurels] We also note that gas from the disrupted host clump 
is initially arranged in a high density spiral, gravitational 
potential of which may well be significant enough to scatter 
planetesimals about and to pump their inclinations. 

This explanation is consistent with the fact that plan- 
etesimals on wider orbits, e.g., a X^ 42 AU, which should be 
affected by the interactions with the gas spiral less, indeed 
have smaller dispersion in their inclination angles. The non- 
zero mean value of i for this population, more distant from 
the planet, should be due to asymmetries developing during 
the disruption process. Comparison with simulation UO be- 
low, which produces far smaller inclinations, demonstrates 
that the asymmetries are probably related to the injection of 
energy in the centre of the host in U15 (and UIO). This en- 
ergy injection increases the entropy of the gas within r < r^j 
around the planet but not at larger radii. The initial en- 
tropy profile is thus strongly unstable to convection in UIO 
and U15, most likely leading to an asymmetric (with respect 
to the initial orbital plane of the clump) expansion of the 
host clump. 

The bottom panel of figure [T] shows the planetesimals in 
the eccentricity - semi-major axis (e-a) plane. The coloured 
curves show the theoretically calculated correlations of e and 
a with slightly different assumptions (explained below) on 
the basis of the approximate theory of §3.4[ All of these 
curves were shifted to ao = 38 AU to fit the simulations 
better. The shift may be empirically justified as following. 
Our approximate theory is based on the assumption of an 
instantaneous disruption of the protoplanet. However, in the 
simulations the protoplanet is not disrupted instantaneously, 
and it does migrate inwards during the disruption, presum- 
ably due to gas gravitational torques. Therefore, in terms 
of our approximate theory, the "effective" location of proto- 
planet disruption should be smaller than the initial radius 
of the host clump's circular orbit. 

Amongst the theoretically computed curves, the green 
dots show the simplest one - equation |22[ We recall that 
this equation was obtained with the help of a Taylor decom- 
position of the specific energy and the specific angular mo- 
mentum in powers of relatively small parameters rjR and 
SvjV . In contrast, the red and the blue curves show the 
predicted orbits of the planetesimals that did not use the 
decomposition. To draw the curves, we consider a ring of 




Figure 7. The inclination (top panel) and eccentricity (bottom 
panel) of the unbound planetesimals versus semi-major axis, a, in 
the simulation U15 at the end of the simulation, t = 10. Individual 
planetesimals' orbital elements are shown with black dots. The 
coloured chevron-shaped curves are theoretical predictions (see 
i ]6.1| . The nearly empty vertical contains many planetesimals still 
bound to the planet, which arc not shown here but analysed in 




planetesimals of radius r = 0.3 and 0.4 AU for the red and 
blue curves, respectively. We also assume that the kick ve- 
locity of the planetesimals after the disruption, Av is equal 
to 5vcirc, where 5 = 0.5 and Vdic is the circular velocity of 
the planetesimal before the disruption (cf. equation^. This 
simple (but somewhat more accurate than the green dots) 
theory also predicts the "V" shaped e — a correlation. 

We also note that our neglect of the second order terms 
in equations [13] and [15] lead to the equation [22] being sym- 
metrical with respect to the sign of a — ao. The more exact 
calculations given with the red and the blue curves are not 
quite symmetrical and also show non-zero eccentricity for all 
the particles. On the other hand, SPH/N-body simulations 
show a wider distribution of eccentricities at a given a, with 
some particles reaching e « 0. Nevertheless, it appears that 
by varying the parameter 5 (which is only constrained to 
be less than unity) and by considering the different rings of 
planetesimals r within the host before its disruption, we can 
qualitatively explain the observed range of e and a as well 
as their correlation. 

The approximate analytical theory of |3.4| also makes a 
prediction on the width of the disruption ring, as u; = 0.5(^, 
where C, is expected to be of order unity. The debris ring in 
simulation U15 has width of ~ 20 AU, commensurate with 
^ = 1. The simple "instantaneous" disruption model is thus 
reasonably accurate in predicting the radial extent of the 
debris disc in this simulation. 
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Figure 8. Same as Figure [7] but at the end of simulation UIO. 



6.2 Unbound debris in UIO 

Figure [S] shows the inchnation and the eccentricity of the 
unbound planetesimals in simulation UIO, along with the- 
oretical curves, in the same format as in fig. W\ Simulation 
UIO differs from U15 by the smaller amount of internal en- 
ergy injected into the host clump by the solid core (planet). 
As the result, the disruption process is much more gradual 
and protracted in UIO than in U15 (see 35|. 

The e-a diagram shows a great deal of similarity be- 
tween simulation UIO and U15, but also some interesting 
differences. For example, the radial width of the disruption 
ring is narrower in the former. This makes intuitive sense 
as there is less internal energy input into the gas clump in 
UIO, and therefore the velocity "kicks" Av after the disrup- 
tion are probably smaller than they are in U15. 

The inclination angle i versus a plot shows somewhat 
larger mean value and larger dispersion of i in UIO than 
in U15. In addition, there are pronounced vertical bands in 
the top panel of Figure IS] Finally, if we break the pattern 
of fig. (tI on the "gas-dominated" central and the "far-out" 
population at a > 43 AU, then its fair to say that in simula- 
tion UIO the far-out population is all turned into the central 
gas-dominated one. 

All of these differences are explainable by the fact that 
the disruption process is more gradual in UIO, which means 
that the tidal tails (arms) of the host clump in the pro- 
cess of disruption (cf. the left panel of fig. jsj survive for 
longer, before being circularised into a diffuse gaseous ring 
(the right panel of fig. [5|. This implies that the planetesi- 
mals are infiuenced stronger by these tidal tails, explaining 
a higher dispersion in inclination i. In addition to that, the 
tidal tails "shepherd" the planetesimals out of the cloud, as 
in seen in the left panels of figures Island [6] causing bunching 
of the orbital parameters in the vertical bands evident in the 
top panel of figure IS] As the spiral arms survive for longer in 



Figure 9. Same as Figs. [7] and Is] except for the red and blue 
curves, as explained in the text (see j ]6.3[ l. Note the much larger 
i Ki population of planetesimals in this simulation compared to 
UIO and U15 



UIO than they do in U15, the bunching effect is far stronger 
in the former than in the latter. 



6.3 Unbound planetesimals in simulation UO 

Figure [9] shows the orbital parameters of the unbound plan- 
etesimals at the end of simulation UO. This figure shows, 
yet again, the familiar e-a correlation diagram but a signif- 
icantly different distribution of planetesimal orbital inclina- 
tions. The red and the blue curves in this case were com- 
puted assuming two rings of planetesimals with r — 0.75 and 
1 AU, respectively, because the planetesimal disc is larger 
in UO than in UIO and U15, and we used S = 0.9 (see j ^6.1[ ). 
The latter is chosen by trial and error to find a visually rea- 
sonable match to the distribution of orbits in the simulation. 

The much stronger bunching of planetesimals towards 
the i = plane is a testament to the different nature of the 
host clump disruption in UO as compared to UIO and U15. 
While in the latter two the inner region of the clump was the 
actual source of the disruption (as the core accretion energy 
was dumped there), in UO the innermost region is "passive". 
Therefore the disruption process, as experienced by the plan- 
etesimals initially located at r ^ 1 AU for this simulation, is 
much less abrupt, leading to even less inclined orbits for the 
unbound particles. In addition, as mentioned in |6.1| the in- 
jection of energy stirs up strong convective motions in those 
simulations, pumping up asymmetries and thus orbital in- 
clinations. 

On the other hand, the larger dispersion in values of i 
in the centre of the planetesimal ring, and the presence of 
the vertical bands in fig. [9] confirms that these features are 
formed by the tidal arms (tails) of the host clump before 
they are wound up and completely erased. 
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r BOUND POPULATION: PLANET 
SATELLITES 

We now switch to discussing planetesimals bound to the 
solid planet, e.g., those that have a negative energy with 
respect to the planet and that are within the planet's Hill's 
radius (see gM. 

Figure [W] presents the orbital parameters of the plan- 
etesimal orbits around the planet for the simulation U15 
(left panel) and UIO (right panel). The red vertical lines 
show the location of the planet's influence radius, r^, de- 
fined by equation |9] This radius divides the region where 
gravitational potential is dominated by the solid planet (in- 
side Ti) and the gas (outside rt). As stated in j |3.3[ the orbits 
of planetesimals are expected to be only mildly perturbed 
within Vi and suffer strong perturbations outside Vi, in fact 
being completely unbound from the planet at r ^ r;. 

Figure [lO] confirms this expectation qualitatively, in- 
cluding the fact that there are very few particles with semi- 
major axis a larger than 0.4 AU. The particles inside r^ tend 
to have mild eccentricity Sj e ^ 0.4, whereas planetesimals 
outside Vi have higher values of e on average. The orbital 
inclination i also increases with a, reaching 10° to 20° out- 
side Ti in simulation U15, and slightly higher values in UIO. 
Note that a larger mean value of i in UIO compared to that 
in U15 found here for the bound population is consistent 
with a similar trend that we found for the unbound popula- 
tion (cf. figures It] and [8| . 

Figure ITT] shows the orbital parameters of the bound 
population of planetesimals in simulation UO. There is a 
marked difference in these distributions compared with that 
for U15 and UIO. Instead of a monotonic increase in eccen- 
tricities and inclinations with increasing a, evident in Figure 
|10[ here there is an almost a step-function change in the na- 
ture of the orbits. Orbits at a ^ 0.3 have small eccentricity, 
^ e ^ 0.2, whereas orbits at a ^ 0.3 have a large spread in 
eccentricities, with some approaching unity. The planetesi- 
mals with the largest values of e may become unbound later 
since their orbits are comparable to r^ ■ There is also a group 
of particles with large inclinations, i ^ 90° at a ~ 0.35 AU. 

The significant differences in the orbits of the bound 
populations of planetesimals between Figures [lO] and [ll] 
demonstrate that the exact way in which the innermost re- 
gion of the host clump is disrupted influences the orbits of 
the "satellites" remaining after the disruption. By extension 
it also means that the results would also be somewhat differ- 
ent if there was a massive gas atmosphere around the solid 
protoplanetary core. 



8 DISCUSSION 

8.1 General conclusions 

In this paper we considered the predictions of the TD hy- 
pothesis for the "planetesimal" debris left over after the dis- 
ruption of a single gas host clump. These gas clumps are the 
parent bodies within which all (or at least most) of planet 
formation action takes place in the TD hypothesis. The most 
important results of our paper can be summarised as follow- 
ing: 

1. In the context of TD hypothesis, large solids form 
only inside the massive gas host clumps (because the den- 
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Figure 11. Same as Fig. |10|but for simulation UO. 



sities of solids initially reflect that of gas, and gas densi- 
ties are orders of magnitude higher inside the host clumps 
than in the "ambient disc"; see, e.g., |Nayakshin|2011a| [Cha| 
& Nayakshin 2011 1 that are then disrupted. The solid de- 



bris is then spread around in a ring around the disruption 
location. TD hypothesis thus generically predicts rings of 
planetesimals remaining after planet formation rather than 
continuous discs. 

2. The role of planetesimals in the formation of planets 
may be very different in the CA and TD scenarios. Whereas 
planets could not have formed without planetesimals form- 
ing first in the CA picture, in the TD hypothesis this is less 
clear. In the latter case, once the initial proto solid core is 
born by gravitational collapse of the grain-dominated inner 
region of the host clump( Nayakshin||2010b 2011al, further 
growth may be dominated by the proto solid core accreting 
~ cm to tens of cm grains. Solids in this size range experi- 
ence aerodynamic drag that is strong enough to dump pos- 
sible centrifugal support and yet small enough to sediment 
to the centre of the host clump within its lifetime (see Fig. 
1). Bodies in the planetesimal size range, i.e., km-sized and 
larger, on the other hand, experience too weak an aerody- 
namic drag. They are not expected to join the solid core en 
masse if they have some orbital support due to an excess 
angular momentum or chaotic convective motions of the gas 
in the centre of the host clump. In this case planetesimals 
are no parent bodies of planets in the strict sense. 

On the other hand, if turbulence, convection or angular 
momentum prevents small grains from reaching the proto 
solid core directly, turning most of these grains into plan- 
etesimals as discussed in §3.2[ then the situation is less clear. 
If this happens sufficiently close to the proto solid core and 
planetesimal densities are high, most of that material can 
be subsequently accreted by the proto core somewhat like 
in the CA theory. In this case planetesimals are also build- 
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Figure 10. The bound population of planetesimals around the solid core (planet). The orbital characteristic arc calculated with respect 
to the planet rather than the star. Left: Simulation U15, Right: simulation UIO. 



ing blocks of planets, with only those on larger less bound 
orbits managing to avoid being accreted onto the core. 

3. Related to the point above, the solid planetary cores 
and the planetesimal debris form at approximately same 
time in our picture. After host gas clump disruption, some 
of the solid debris around the planetary core could fall on 
it and be accreted in direct collisions, so that planetary 
growth could actually continue in a manner analogous to 
that of the later oligarchic stages of solid core assembly in 
the traditional scenario (e.g., Safronov|1972 I. In particular, 
as velocity dispersion of the disrupted population is high, 
e.g. a fraction of a km s~^, only large planetesimals of size 

^ 1000 km would continue their growth, whereas smaller 
objects would accrete on large bodies or fragment in colli- 
sions with smaller bodies. If they fragment to sizes as small 
as a few metres then their inward radial migration inside the 
gas disc becomes important. 

4. Also related to point (2) above, the mass budget of 
planetesimals can be potentially very different in the two 
planet formation theories. In the CA scenario, the initial 
population of planetesimals must have had a total mass 
larger than the total high-Z element (other than H and He) 
mass of the resulting planets. This does not have to be the 
case in the TD scenario. In particular, there does not seem 
anything wrong with assuming that the debris population 
may be far less massive than the protoplanetary core. 

5. In the CA scenario, the properties of the protoplane- 
tary disc are expected to be a rather smoothly varying func- 
tion of radius, except for the regions near the ice line (e.g., 
Armitage|2010 1. Thus one expects that planetesimal proper- 
ties are also a smoothly varying function of R. In contrast, in 
the TD hypothesis, it is the host clump that determines the 
outcome of the planet formation process. The internal struc- 
ture of the clump is a strong function of its mass ( |Nayak-| 
shin 2010b 2011a I; that structure is non-uniquely coupled 



to radius 7? where the clump is disrupted. It is possible 
that clumps of different history, mass, internal structure are 
disrupted at roughly same location. Further, we envisage a 
number of host clumps forming and being destroyed during 
the early gas-rich phase as required to explain the FU Ori 
outbursts of young protostars ( [Nayakshin fc Lodato|[2011[ ). 
This implies a far greater diversity in the properties of the 
solid debris at the same spatial location than is possible in 
the CA model. 

6. We found in this paper that disruption of a host 
gas clump generically produces a "V" -shaped pattern in the 
eccentricity-semi-major axis space for the planetesimal pop- 
ulation released by the disruption. The planetesimals in the 
centre of the pattern have nearly circular orbits while whose 
at the edges of it are much more eccentric. The maximum 
eccentricities possible can be estimated based on the simple 
analytical theory; from equations [22] to |25| it follows that, 
within a factor of ~ 2, it is emax ~ rn/R ~ 0.2 for host 
clump mass of ~ 10 Jupiter masses. 

7. The eccentricities and inclinations of the disrupted 
planetesimal population are generally much too high to al- 
low planetesimals to stick to one another. The typical disper- 
sion velocity of the disrupted population is 5v ~ vkTh/ R ~ 
Q.Ivk- At the distance of the Kuiper belt, for example, this 
yields dispersion velocity of the order of 300 m/s. Only 
Pluto-sized bodies could survive equal-size collisions at these 
velocities. Accretion on massive solid cores, however, is al- 
lowed as already noted in (3) above, so the post-disruption 
evolution may be somewhat similar to the final phases of 
the run-away growth of planetary embryos in the traditional 
scenario 



(e.g. 



Safronov 1972}. 



8. We also find that solid bodies on orbits tightly bound 
to the planetary core could survive the disruption of the 
gas host clump, remaining bound to the planet. These bod- 
ies are satellites of the planetary cores. Satellites on orbits 
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more strongly bound to the planet are affected by the gas 
envelope destruction less than those on less bound orbits. 
Some of the outermost satellites may even find themselves 
in eccentric counter-rotating orbits (e.g., see fig. fOl. This 
is qualitatively similar to the satellites of the giant planets 
in the Solar System, although we note that our simulations 
are not designed to study these issues; one also would need 
to model composition, size differences and long-term surviv- 
ability of the bound objects. 



8.2 Potential relevance to the Solar System 

Even though we did not specifically attempt to reproduce 
the structure of the outer Solar System here (which would 
require an additional study of how the system evolves on 
the 4.5 Giga year time scale), our simulations do have po- 
tentially interesting implications for it. We mention these 
implications only briefly here, with the intent of quantifying 
them in a future publication: 

• One of the key properties of the trans-Neptunian region 
of the Solar System, including the Kuiper belt, is the "mass 
deficit problem" . The current mass of the Kuiper belt is es- 
timated at ~ 0.01 - 0.1 M® ( [Bernstein et al.|2004[ ). On the 
other hand, to grow the observed populations of solid bodies 



in the context of the Safronov ( 1972 1 model for formation 
of solids, 10 to 100 M^ of solids is required (e.g., [Kenyon| 
' fc Luu|l999"| . Thus, removal of > 99.9% of solid material is 
required. In the TD hypothesis, solids are made within the 
central fractions of AU of the host gas clumps ( |Nayakshin| 
2010b I. The host clumps contain as much as ~ 60 M® of 
high-Z elements ( Nayakshin|2011a | which are then partially 
used to make the massive solid cores, the "planetesimals" , 
and partially remain bound to gas in small grains. If it is 
possible to build solid cores as massive as ~ 10 M^ by di- 
rect gravitational collapse and accretion of ~cm size grains 
( Nayakshin|2011a I , it should also be possible to build Pluto- 
sized objects without having to put tens of M^ of material 
into the planetesimals. There is thus no mass deficit problem 
for the Kuiper belt in the TD scenario. 

• There is a sharp outer edge to the Kuiper belt at around 
7? ~ 50 AU which is currently not well understood (e.g., 
Morbidelli et al.]|2008 1 . The disruption of a gas host clump 
naturally produces a ring with sharp outer and inner edges 
(point 1 above) because the range of orbits available to plan- 
etesimals after disruption is limited by the parameters of the 
initial host clump. This could naturally account for the outer 
edge of the Kuiper belt objects. 

• The structure of the classical Kuiper belt is best ex- 
plained by our model if we assume that the left hand side 
(lower a) of the disruption ring (e.g.. Fig. S| has been 
destroyed by interactions with planets. Realistically, the 
planet, left at a ~ 35 — 39 AU at the end of our simulations, 
could have continued to migrate inward via type I migration 
in the gas disc: the disrupted gas ring is still quite massive 
(5 Jupiter masses for our simulations here, and there would 
be more at larger R to effect the migration of the host clump 
in the first place). Using the standard type-I gas migration 



fast to allow Neptune to migrate inward significantly to end 
up, for example, in the configuration proposed by the NICE 
model of the outer Solar System ( Gomes et al.|2005 l. While 
migrating inward in the gas disc, Neptune would have scat- 
tered the left part of the "V" -pattern of the planetesimals, 
but leaving behind the right hand side of the pattern as a 
belt reminiscent of the Kuiper belt. 

• Kuiper belt contains the hot and the cold populations. 
These are best accounted for by two different gas clump 
disruptions in our model. The hot population in our model 
would have most naturally resulted if the host gas proto- 
planet rotation axis were highly inclined to the orbital plane. 
In this case one could end up with higher inclinations than 
we obtained here. 



Nesvorny et al. (20101 show that if planetesimals form 



by a local gravitational collapse in a high density environ- 
ment, then one can naturally explain the surprisingly large 
fraction of binaries in the ~100 km class low-inclination ob- 
jects in the classical Kuiper Belt. We note that our model 
produces "planetesimals" in a similar fashion (gravitational 
collapse in a high density environment) albeit inside of the 
host clumps rather than the disc. Therefore by extension 
one might expect to see a high fraction of massive objects 
locked in binaries after the dissipation of the host clump. 
This idea however needs to be checked with a longer time 
A'^-body calculation as there are non trivial physical con- 
straints on coUisional destruction of binaries in the Kuiper 
Belt ( jNesvorny et al.||2011[ ). 



9 CONCLUSIONS 

We have considered the origin of solid debris, such as comets, 
asteroids and large bodies such as Pluto in the context of a 
recent planet formation hypothesis (Tidal Downsizing). We 
assumed (see \ 3.2 1 that the solid "debris" is formed in a way 
very similar to the massive solid protoplanetary cores them- 
selves, e.g., inside the ~ 10 Jupiter mass host gas clouds. 
The latter are eventually disrupted either tidally or due to 
an internal energy release during the solid core accretion. 
The release of these solid bodies into the field forms rings 
potentially reminiscent of the Kuiper belt and the debris 
discs around nearby main sequence stars. While much work 
remains to be done to detail predictions of the TD hypoth- 
esis further, it is already clear that these predictions are 
sufficiently different from the standard planetesimal-based 
paradigm for planet formation ( Safronovl|1972 | to be criti- 
cally tested by observations in the near future. 

As an astrophysical aside, we note that the rapid accre- 
tion of large solid bodies (Pluto-like and even Neptune-like) 
in the TD scheme suggests that planet and solid debris for- 
mation is a very robust process and may even occur in very 
crowded environments such as the inner parsecs of galax- 
ies (Nayakshi n et aL||2012 Zubovas et al.pOll I; this is an 
untenable proposition in the CA model. 



(e.g., Tanaka et al. 20021 for Neptune initially located at 
« 40 AU, assuming the disc mass being ~ 10 Jupiter masses, 
and the gas disc aspect ratio of H/R ~ 0.1 yields type I 
migration rate time scale of ~ 10^ yrs. This is sufficiently 
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